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Abstract 

We compare static arbitrage price bounds on basket calls, i.e. bounds that only 
involve buy-and-hold trading strategies, with the price range obtained within a multi- 
W | variate generalization of the |BS73| model. While there is no gap between these two 

^«-^ ' sets of prices in the univariate case, we observe here that contrary to our intuition 

about model risk for at-the-money calls, there is a somewhat large gap between model 
prices and static arbitrage prices, hence a similarly large set of prices on which a mul- 
tivariate }BS73| model cannot be calibrated but where no conclusion can be drawn on 
■ the presence or not of a static arbitrage opportunity. 
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In the classic unidimensional |BS73j framework, there is no gap between the range of model 
prices on one hand and the range of prices that create a buy-and-hold arbitrage opportunity 
on the other. This means in practice that if we can't extract a market implied volatility from 
the price of a call, then we know that there either an error in the data or a buy-and-hold 
arbitrage opportunity in the market that makes this price unviable on the long run. 

The [BS73] call price formula is a strictly increasing function of the volatility and can 
be easily inverted, we also know from |LL00j for example that testing for the presence of a 
buy-and-hold arbitrage is equivalent to testing for positivity, monotonicity and convexity on 
call prices, which can be done by inspection. The unidimensional setting is then extremely 
favorable, since model calibration can be achieved at very little numerical cost, and a final 
answer is obtained as easily on the presence or not of buy-and-hold (or static) arbitrage 
opportunities. Because of these numerical properties and for consistency, call prices are 
always quoted in terms of their implied volatility. 

Unfortunately, those key numerical properties are lost in a multivariate setting. As we 
will see below, calibrating a multivariate [BS73 j model is a non trivial exercise and testing 

*EECS Dept., U.C. Berkeley, CA 94720. alexandre.daspremont@m4x.org 



1 



for the absence of static arbitrage between basket options becomes an NP-Hard problem (see 

[EEl). 

Equity derivatives markets don't make very frequent use of basket options beyond some 
elementary spread options. However, interest rate derivatives market volatility information 
is mostly concentrated in caps and swaptions, which can be seen as basket options on for- 
ward rates (see [Rcb98 or jd'A03c| among others). While empirical evidence suggests that 
multivariate lognormal approximations to market models (see [BDB99] or |d'A03cj ) calibrate 
very well to prices of caps and swaptions, a satisfactory joint model of correlation and smile 
features has yet to be designed. This means in practice that there remains a range of option 
prices on which these models cannot be calibrated but where no conclusion can be drawn on 
the presence or not of a static arbitrage. In this work, we try to quantify the magnitude of 
this range of prices. 

The paper is organized as follows, in section two we show how to derive bounds on the 
price of a basket call within a multivariate lognormal model. In section three we detail 
various relaxation techniques to compute static arbitrage bounds on the price of baskets 
while section four details some numerical results. 



2 Model price bounds 

In this section, we compute the range of prices covered by a multivariate lognormal model 
when calibrated to a set of (liquid) market instruments. This is a hard numerical problem 
in general but excellent estimates of these bounds can be computed in a multivariate [BS73 
model. We briefly describe the results of |d'A03c| in a simple equity setting, for more details 
on market model approximations and interest rate options pricing, we refer the reader to 
|Reb98| . |BDB99| or |d'A03c| . 

In this setting, the dynamics of the assets F\ are given by: 

dF\ = FydW s , 

where W s is a n-dimensional Brownian motion and a % G R n for % = 1 , . . . , n are the volatility 
parameters. We shall denote by X e S n the corresponding covariance matrix, with = 
a %T (ji . The sum of lognormally distributed assets is not lognormal, but we sue a lognormal 
approximation to price basket calls. From |d' A03cj . we know that the price of a basket call 
with payoff: 



WiFl - K 




at time T, can be approximated by a |BS73j call price using an appropriate variance Vt such 
that: 

C = BS(w T F u K, T, V T ) = (w T F t )Af(h(V T )) - KM (h{V T ) - ^V^ , (1) 



2 



where M{x) is the CDF of the normal distribution: 



(in(^) + iy, 



with 

Vr = Tr (SIX) T, 
where f2 e S n is the matrix given by SI = ww T and 

Wi - 



w 



T Ff ; 



This gives swaptions prices accurate to within 1-4 basis points in the Libor market model. 
Furthermore, a good estimate of the hedging tracking error can be computed by robustness 
argument (see |d'A03b| ). 



Since the |BS73j formula is strictly increasing with its variance term Vt, computing 
bounds on the model price of a basket call given market price data on other baskets (with 
the same maturity) is equivalent to solving the following semidefinite program: 

max./min. Tr(f2 X) 

subject to Tr(fiiX) = V T) i, i = l,...,m (2) 

x y o, 

in the variable X e S n with Vr t i such that: 

BS(wfF t , Ki, T, V ijT ) =pi, i = 1, . . . , m, 

where pi are the market prices of basket call options with weights Wi, maturity T and strike 
Ki. Note that we implicitly assume that all the options have the same maturity, and that, 
without loss of generality, the risk-free interest rate is zero (we compare prices in the forward 
market). Let us remark however that the multiperiod generalization has exactly the same 
format (see |d'A03c| ). 

This last semidefinite program can be solved very efficiently using algorithms such as the 
one by [Stu99j, we refer the reader to |NN94j or |BV04j for further details. This means that 



given a certain number of prices of liquid market instruments, we can efficiently compute 
upper and lower bounds on the model price of another instrument. 

If the market price of an instrument falls outside of these bounds, we then know that 
calibrating the model to this instrument together with the original set of prices will be 
infeasible. This means that there is a dynamic arbitrage opportunity if (and only if) the 
market dynamics follow that of the model. In practice however, this more often means either 
that the model dynamics are not rich enough to capture the market price features or that 
there is a problem with the market data set (liquidity, outliers, missing data, etc). In the 
next section, we discuss ways of refining this diagnostic. 
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3 Static arbitrage bounds 



In the previous section, we computed bounds on basket option prices assuming that the 
underlying assets followed lognormal dynamics. In this section, we are looking for price 
bounds on options (given prices of other liquid options) without any assumption on the 
asset distribution. The range of prices covered will of course be much larger than in the 
lognormal case but option prices falling outside of this range generate buy-and-hold arbitrage 
opportunities, which are much more robust to liquidity issues than the model arbitrage 
detailed in the last section. 

Let p e R+ , K e R™, w e R n , wi e R n , i = 1, . . . , m and K > 0. We consider here 
the problem of computing upper and lower bounds on the price of an European basket call 
option with strike K and weight vector Wq: 

min./max. E 7r (w^x — K ) + , . 

subject to ~E n (wfx — Ki) + — p i: % = 1, . . . , m, 

with respect to all probability distributions 7r with support in R™ on the asset price vector 
x, consistent with a given set of observed prices Pi of options on other baskets. Note that we 
again implicitly assume that all the options have the same maturity, and that, without loss 
of generality, the risk-free interest rate is zero (we compare prices in the forward market). 

We know from BP02J that the problem described in (j3J) is NP-Hard and in the next 
section we describe several relaxation techniques providing upper (resp. lower) bounds on 
the upper (resp. lower) bounds described by (J3J). 

3.1 Linear programming relaxation 

In dimension one we know that if C(K) is a function giving the price of an option of strike 
K, then C(K) must be positive, decreasing and convex. With C(0) = S, we have a set of 
necessary conditions on call prices for the absence of arbitrage. 

In fact it is well known (see [LLQOj or |BP02j among others) that these conditions are 
also sufficient, so there is no arbitrage between some given market call prices Pi if and only 
if there is a function C(K) such that: 

• C(K) positive 

• C(K) decreasing 

• C(K) convex 

• C(Ki) = pi and C(0) = S, i — 1, . . . , m. 

These conditions are easily generalized to the multidimensional case as follows: given a set 
of market prices for basket calls C(w{, Kj) = pi and suppose there is no arbitrage, then the 
function C(w,K) must satisfy: 

• C(w,K) positive 
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• C(w, K) decreasing in K, increasing in w 

• C(w, K) jointly convex in (w, K) 

• C(wi, Ki) = pi and C(0) = S, i = 1, . . . , m. 

The key difference is here that these conditions are not sufficient. Nevertheless, we can form 
a relaxation of problem ©as: 

min./max. C (wq, K ) 

subject to C(w, K) positive 

C(w,K) decreasing in K, increasing in w (4) 

C(w, K) jointly convex in (w, K) 

C(wi, Ki) = pi and C(w, 0) = w T S, i — 1, . . . , m. 

In this form, this is an infinite dimensional linear program and not directly tractable. How- 
ever, we can discretize (j3J into a finite linear program by sampling the constraints on the 
data points (wi,Ki). We get the following program 

maximize / minimize p$ 

subject to gf((wj, Kj) - {w h Ki)) < pj - p h i,j = 0,...,m ,^ 

9i,j > 0, -1 < gi, n +i < 0, i = 0, ...,m, j = l,...,n 
gf((wi,Ki)) = p i} i = Q,...,m, 

in the variables po E R + and gi G R n+1 for i — 0, ...,m. This is a linear program with 
m x (n + 1) variables and can be solved very efficiently Furthermore, [dEG02j show that 
program (j3J) is in fact an exact discretization of program (@J) and that an optimal solution 
to (jlj) can be constructed from that of 

This first relaxation technique allows us to get upper and lower bounds on the solution 
to at a minimal numerical cost, clEG02j show that these bounds are exact in some 
particular cases, but in general nothing can be guaranteed about their performance. 

Necessary and sufficient conditions describing when a function can be written as the price 

p = 'E 7T (w T x - K) + 

of a basket call with weight w and strike K were derived by [HS90 in their investigation of 
production functions. They get the following result: a function can be written 

C(w, K)= [ (w T x - K) + dn(x), 
Jr.™ 

with w G R™ and K > 0, if and only if: 

• C(w,K) is convex and homogenous of degree one 

• for every w G R™ + , we have lim^^oo C(w, K) = and lim^^ + dC j^ K ^ = — 1 
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. F(w) = f °° e~ K d belongs to C^(n%) 

• For some w G R™, (— l) k+1 D^...D^ k F(Xw) > 0, for all positive integers k and A G R++ 
and all £i> • • • > £fc in R™ . 

The key difficulty here is that the last two smoothness and monotonicity conditions are 
numerically intractable, so this result is of very little practical help in refining the conditions 
used in (j5J). 

3.2 A moment approach 

In this section, we look for ways of improving the relaxation technique in (|SJ). NP hardness 
means that we can't hope to get a tractable set of necessary conditions to solve the problem 
exactly, here instead we look for additional conditions on prices that produce a sequence of 
successively tighter bounds on the solution to (j3J). 

The integral transform approach above suggests a link to moment theory. In fact, as 
detailed in |d'A03aj . numerically tractable conditions for the existence of a measure 7r such 
that p = E 7T (w T x — K) + can be obtained by a generalization of Bernstein-Bochner type 
results to the payoff semigroup (see [BCR84 for a complete exposition). We briefly recall 
this construction below. 

We suppose that the market is composed of cash and n underlying assets Xi for i = 
1, . . . , n with x G R™ . We suppose that the forward prices of the assets are known and given 
by pi, for i — 1, . . . , n, hence Wi is the Euclidean basis and K j = for i — 1, . . . , n. In addition 
to these basic products, there are m + 1 basket straddles on the assets x, with payoff given 
by \w^ +i x — K n+i \, % = 1, . . . , m. Because a straddle is obtained as the sum of a call and a 
put, we get the market price of straddles from those of basket calls and forward contracts by 
call-put parity and the static arbitrage problem on straddles is strictly equivalent to problem 
(J2J) since we always assume that forward prices are quoted in the market. 

For simplicity, we will note these payoff functions e^, for % = 0, . . . , m + n, with cq{x) = 
\wqX — K \, ei(x) — Xi for i = 1, . . . , n and e( n+ j)(x) = \wfx — for j = 1, . . . , m. In what 
follows, we will focus on the commutative semigroup (§, •) generated by the payoffs ei(x) for 
i = 0, . . . , m + n, the cash 1§ and their products. 

S = {l,xi, . . . , Iw^x - K m \,xl, . . .,Xi\wjx - Kj\, . . .} (6) 

Let us recall that a function / : S — > R is called positive semidefinite iff for all finite 
families {sj} of elements of §, the matrix with coefficients f(siSj) is positive semidefinite. 
We then get the following result from [d'A03a , suppose the asset distribution has compact 
support K and S is the payoff semigroup defined above. A function f(s) : § — > R can be 
represented as 

f{s) =E„[s(z)], forallsGS, (7) 
for some measure v on K, and satisfies the price constraints in © if and only if: 

(i) f(s) is positive semidefinite, 
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(ii) /(ejs) is positive semidefinite for i = 0, . . . , n + m, 

(iii) (f3f(s) — Yli^o 1 f( e i s )) 1S positive semidefinite, 

(iv) f(ei) = Pi for % = 1, . . . , n + m. 

Furthermore, for each function / satisfying conditions (i) to (iv), the measure v. A similar 
result holds in the case where the support of v is not compact (see |d'A03aj ). 

The above result shows that testing for the absence of static arbitrage between the 
securities in S, i.e. the set of straddles and their products, is equivalent to testing the 
positivity of an infinite number of matrices. This gives a direct recipe for writing a relaxation 
to problem (jHJ): by considering a finite subset of these matrix inequalities, © can be turned 
into a semidefinite program. We summarize this procedure below. 

We begin by recalling the construction of moment matrices in jLasOlj . Let A(S) denote 
the real algebra generated by the payoffs in §. We adopt the following multiindex notation 
for monomials in *4(§): 

e (x) := e ° (x)e l 1 (x) ■ ■ ■ c m + n (x), 

and we let 

Ve = (1) GOj • • • > Cm+ri) e e l; • • • j e 0' * * * J e m+n) VV 

be the vector of all monomials in ^4(S), up to degree d, listed in graded lexicographic order. 
We note s(d) the size of the vector y e . Let y G R S ( M ) be the vector of moments (indexed as 
in y e ) of some probability measure v with support in R", we note M d (y) G R s ( d ) xs ( d ) ) the 
symmetric matrix: 

M d (y)ij= (y e ) i (x){y e ) j {x)diy(x), for i,j = 1, s{d). 

In the rest of the paper, we will always implicitly assume that y± = 1. With /3(i) the exponent 
of the monomial (y x )- and conversely, i(/3) the index of the monomial e 13 in y e . We notice 
that for a given moment vector y G Hf^ ordered as in the first row and columns of the 
matrix M d (y) are then equal to y. The rest of the matrix is constructed according to: 

M d (y)ij = yua+p) if M d (y) iA = y i{a) and M d (y) hj = y m . 

Similarly, let g G *A(§), we derive the moment matrix for the measure g(x)dv on R", 
noted M d (gy) G S s ^ d \ from the matrix of moments M d (y) by: 

M d (gy) i;j = / (y e ) i (x) (y e ) j (x)g{x)dv{x) for i,j = 1, s(d). 

JR™ 

The coefficients of the matrix M m (gy) are then given by: 

M d {gy) itJ = ^2g a yi((3(i )+ (3( j)+a ) (9) 

Q 
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We can now form a semidefinite program to compute a lower bound on the optimal 
solution to (J3J) using a subset of the moment constraints in (J7J), taking only monomials and 
moments in y up to a certain degree. Let N be a positive integer and y G ~R?( 2N \ a lower 
bound on the optimal value of: 

minimize p := E^[e (a;)] 

subject to E v [ej(x)] = Pi, i — 1, . . . , n + m, 

can be computed as the solution of the following semidefinite program (again see |NN94j or 

IHYOU): 

minimize y<± 
subject to M N (y) >z 

M N {e jy )hO, for j = 1, . . . , n, (10) 

y(j+2) — Pj, for j = 1, . . . , n + m and s G § 

where s is such that i(s) < s(2N). As N increases, program (jl()j) provides an increasingly 
precise relaxation of problem ©. 



4 Numerical results 

Here, we try to quantify on simple examples the magnitude of the gap between model prices, 
the price range obtained using the various relaxation techniques detailed above and the exact 
price bounds in Q. 

4.1 Discrete model 

Here, we simulate a set of arbitrage free basket call prices using a simple discrete model. 
Given these prices and the absence of arbitrage between basket calls, we study the price 
bounds induced on another basket call. 

Of course we can't compare these bounds with the exact solution, however we can try to 
compare these bounds with inner bounds obtained by maximizing and minimizing the price 
C(cuo, Kq) over a set of probability measures satisfying the price constraints C(uii, Ki) = p^. 
If we consider only discrete measures, this becomes a (very large) linear program. 

We suppose here that the asset price at maturity T lies within the unit box [0, 1]™. We 
then discretize the probability density using a grid with iV bins per asset. The problem of 
finding (inner) upper and lower bounds on a basket (u , K ) can be written as: 

max./min. E u (u%x-K ) + 

subject to Yj u {ujJ x — K^j = pi, i = l,...,m, 

which is a linear program of (exponential) size N n in the discrete measure v. 
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We test these upper and lower bounds in dimension two. Program can be written: 

max./min. Y,k,i=o,...,N MW-^, l/N] T - K ) + 

subject to Y J k,i=o,...,N u ki(u?[k/N,l/N] T - Ki) + = pi, i = l,...,m. 

which is a linear program in the variable v G R iVxAr . The assets are noted x\, x% and we look 
for bounds on the price of an index option with payoff {x% + x 2 — K) + . To produce price 
data, we use a simple discrete model for the assets, their distribution has finite support and 
is given by: 

x = {(0, 0), (0, .8), (.8, .3), (.6, .6), (.1, .4), (1, 1)} 

with probability 

{.2,.2,.2,.1,.1,.2}. 

The input data set is composed of the forward prices together with the following call prices: 

(.2xi + x 2 - .1)+, (.5a;i + -8x 2 - .8)+, (.5xi + -3x 2 - .4)+, 

(xi + .3x 2 - .5)+, (xi + .5x 2 - .5)+, (xi + .4x 2 - 1)+, (x 1 + .6x 2 - 1.2) + . 

We plot the inner and outer bounds obtained using this data in figure (JTJ). We observe that 
sometimes the bounds match, i.e. the price bounds given by the relaxation are tight, while 
sometime there is a gap and not much can be said about the relaxation's suboptimality. 

We also examine how these bounds evolve as more and more instruments are incorporated 
into the data set. For a particular choice of strike price (here K = 1), we compute the outer 
bounds © and inner bounds obtained when using only the k first instruments in the 
data set, for k = 2, . . . , 7. The result is plotted in figure (J2J). We notice that it takes relatively 
few basket prices to get good bounds on the price of the target option. 



4.2 Multivariate lognormal model 



Here again, we simulate a set of arbitrage free basket call prices using this time the multi- 
variate jBS731 model detailed in section one. Given these prices and the absence of arbitrage 
between basket calls, we study the price bounds induced on another basket call. 

We look for bounds on the price of an index option: (.2 Yui=i x * — ^)+> given the price of 
all at-the-money single asset calls and at-the-money basket calls with the following weights: 

0.33 0.33 0.33 0.00 0.00 
0.00 0.00 0.33 0.33 0.33 
0.40 0.20 0.20 0.20 0.00 

The assets initial values F\ are (0.03, 0.03, .05, .07, .07) and the model covariance matrix 
X is given by: 

/ 0.06 0.04 0.04 0.04 0.04 \ 

0.04 0.06 0.04 0.04 0.04 

0.04 0.04 0.06 0.04 0.04 

0.04 0.04 0.04 0.06 0.04 

\ 0.04 0.04 0.04 0.04 0.06 / 
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Price bounds 
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0.4 0.6 



strike 

Figure 1: Comparison of the inner bounds computed by discretization (solid lines, computed 
using (JUJ) and the outer bounds obtained by relaxation (dashed lines, computed using ©). 
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number of instruments 

Figure 2: Inner bounds (solid lines, computed using (JUJ)) and outer bounds (dotted lines, 
computed using (|SJ)) versus number of instruments in the data set. 
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We plot the inner and outer bounds obtained using this data in figure (j3J). Since the 
[BS73] of an at-the-money call is not very sensitive to the volatility, we could have expected 
the call price near the money to be somewhat insensitive to model specification. The fact 
that the gap between the inner model price bounds and the outer static arbitrage bounds is 
so large seems to directly contradict this intuition, showing that in fact the range of arbitrage 
free prices for basket calls is much larger than the range of prices that can be attained by a 
multivariate |BS73| model. 

4.3 Moment constraints 

For numerical reasons, we consider a model with two assets x%,X2 and look for bounds on 
the price of the basket \x\ + x 2 — K\. We use a simple discrete model for the assets: 

x = {(0,0), (0,3), (3,0), (1,2), (5, 4)} 

, with probability 

p= (.2, .2, .2, .3,1) 
, to simulate market prices for the forwards and the following straddles: 

\xi — .9|, \xi — 1|, \x2 — 1.9|, \x2 — 2|, \x2 — 2. 1 1 . 

The results are detailed in figure Unfortunately the code from [Stu99j, although very 
robust in general, is not as stable on the moment problems detailed here. A different im- 
plementation using large scale optimization techniques such as spectral bundle should prove 
more appropriate here. 

5 Conclusion 

We have detailed some tractable relaxation procedures to test for the absence of static arbi- 
trage between basket calls. A comparison with the range of prices attained by a multivariate 
[BS73j model shows that this simple model only covers a relatively small portion of the 
range of arbitrage free prices. Running the same kind of test on a model with a richer 
smile structure, if numerically feasible, would be very interesting. Also, the numerical issues 
encountered in the moment relaxation technique are very surprising and probably deserve 
more investigation. 
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Figure 3: Comparison of the inner bounds computed using the Black & Scholes model 
(solid lines, computed using (J2J) and the outer bounds obtained by relaxation (dashed lines, 
computed using (j3J)). 
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Figure 4: Upper and lower price bounds on a straddle (solid lines). The dashed lines represent 
the payoff function and the actual price. 
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